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I. INTRODUCTION 

In a renormalizable quantum field theory at 4- dimensional Euclidean space time, the 
local interaction of a fermion system with external fields is bilinear in fermion fields and 
realized in terms of Dirac operator D, 

D = y - s + ipj 5 = <9 M - iv^ - ia M 7 5 = -V], , (1) 

where s, p, v^, a M are hermitian fields fermion coupled with and 7 matrix is hermitian. We 
have decomposed these fields in terms of their Lorentz structures 0. In terms of Dirac spinor 
i/j and ip, the action of renormalizable interaction is 



* Mailing address 

1 In general, there should be tensor field, but for simplicity of discussion, in this paper, we donot 
consider the tensor field. 
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J d 4 x ipDip 

This action is invariant under following local chiral transformations 



(2) 



iP(x) -> iP'(x) = [V L {x)P L + V R (x)P R ]iP(x) 

J(x) - J\x) = [V L (x)P R + V R {x)P L ][J{x) + x ][Vl(x)P L + V R (x)P R ] (3) 

with Vl{x) and V R (x) be left and right chiral rotation matrices and project operator Pr = 
|(1 ± 7s). The external field J(x) is defined as 

J(x) = —if [x) — ijt, (x)7 5 — s[x) + ip(x)j 5 (4) 

Suppose we have Nf Dirac spinors in our fermion system, this local chiral symmetry then 
is U L (N f )®U R (N f ). 

The contribution of action (Q) to a path integral of 4- dimension space time quantum 
field theory can be written as 

where I and I are external sources of the fermion system. D is Dirac operaor which in 
general depend on external fields. Due to bilinear property of fermion interaction, we 
have integrated out fermions in the path integral ([|) and found that the contribution of 
fermion to the path integral is realized through fermion determinant TrlnD and propagator 
D~ l (x,y) = (xlD^ly) in presence of external fields. Once these two objects are obtained, 
the contribution of fermion system to a 4-dim space time renormalizable field theory is 
known. 

In practice, for the physical needs and regularization of infrared divergence, a constant 
bare mass m is usually added into the theory which is equivalent to extract out from the 
scalar field s a condensation part s — > s — m. With this quark mass term, interaction action 
(0) become 

J d A x tjj(D + m)i) (6) 

and two ingredients of fermion system in ([|) become logarithm of fermion determinant 
Trln(D + m) and fermion propagator (D + m)~ 1 (x, y). 

Adding in quark mass term explicitly violate original chiral symmetry from Uz,(Nf) (g) 
U R (Nf) to Uv{Nf), whose transformation is in fact requiring the original left handed and 
right handed rotations be same Vl(x) = V R (x) = h\x). 
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With constant quark mass term, a way to keep the theory still be chiral symmetric is 
nonlinear realization of the theory, i.e., we need to introduce into theory a local field Q(x) 
which transform nonlinearly under local chiral symmetry (|3|) as 

Q{x) -> Q'(x) = h\x)VL(x)Vl(x) = V R (x)tt(x)h(x) (7) 

and rotated external fields 

J n (x) = [Q(x)P R + rt(x)P L ] [J(x) + ip\ [Q(x)P R + Q\x)P L ] 

= ~sn{x) + ipnix)^ + j> n (x) + </l q(x)7 5 (8) 

in which h(x) is a transformation matrix for induced hidden local symmetry Uy(Nf). Re- 
placing the external fields in D + m in (Q) with Dq + m, i.e. 

D + m ^ D n + m (9) 

with Dq = y q — sq + ipn^, = (9 M — iv^ — ia^s- One can easily check that under local 
chiral symmetry Ui{Nf) <g> U R (Nf) transformation (^|) and (0), 

J n (x) ^ = /i+(x)[J n (x) +i i 9 a ,]/i(x) (10) 

D n + m->D| J + m = /i t (a:)[Dn + m]/i(a;) (11) 

i.e., they are covariant quantities. In terms of rotated quantities, local chiral symmetry 
Uz,(Nf) ® U R (Nf) is nonlinearly realized through hidden symmetry Uv{Nf) by /i(a;). Once 
our formulation keep this local hidden symmetry Uy(Nf), it automatically keep original local 
chiral symmetry Ul{Nj) <g> U R {Nf). If our formulation keep the local symmetry Uy(Nf), it 
then keep local chiral symmetry Ul{Nj) ® U R (Nf) in terms of f2 field. So with constant 
quark mass term, to keep local chiral symmetry UL{Nf)®U R (Nf) in the theory, we only need 
to replace original sources with the rotated sources. In this sense, one can easily generalize 
the theory with symmetry Uy(Nf) to the case of symmetry Ul(Nj) eg) U R (Nf) in terms of Q 
field. We in this paper will limit ourselves with symmetry Uy(Nf) and unrotated sources. 

Beyond constant quark mass, in the literature, there is discussion on how to generalize the 
constant fermion mass to a matrix jlj] to respect further breaking of the Uy(Nf) symmetry in 
the fermion system. In this work, we are interested in the generalization in another direction- 
the case of dynamical violating the chiral symmetry. Within context of dynamical chiral 
symmetry breaking (DCSB), the most general situation is not the appearance of a hard 
fermion mass m, but a momentum k 2 dependent fermion self energy S(fc 2 ). In fact, in the 
quantum field theory, even the chiral symmetry is broken by a momentum independent hard 
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mass term at leading order of calculation, include in high order quantum corrections, the hard 
mass will be replaced by momentum dependent fermion self energy. Taking hard fermion 
mass often cause extra ultraviolet divergence, instead, momentum dependent fermion self 
energy usually will suppress the ultraviolet divergence. For example, in the case of massless 
QCD, the fermion self energy damping at least as 1/k 2 at ultraviolet momentum region, if we 
take a hard fermion mass to substitute self energy, we will over estimate its contribution to 
physics and cause extra divergence which usually lead many physical fine tuning problems. 
To respect the momentum dependence of fermion self energy in the theory, in this paper, we 
generalize the conventional Schwinger proper time formulation to deal with corresponding 
fermion determinant and propagator. For simplicity and as the first step of research in this 
direction we only discuss the real part of fermion determinant known conventionally as 
effective action and part of coincidence limit of fermion propagator related to fermion pair 
condensation. The imaginary part of fermion determinant related to anomaly and remaining 
nonlocal part of fermion propagator will be discussed else where. 

This paper is organized as follows: In section II, we review the conventional Schwinger 
proper time formulation for the real part of fermion determinant and coincidence limit of 
fermion propagator. In section III, we discuss how to generalize the formulation from a hard 
mass case to momentum dependent dynamical fermion self energy. Section IV is responsible 
for the calculation of the real part of fermion determinant and section V is responsible for 
the calculation of the coincidence limit of fermion propagator in which the tensor structure 
is ignored in our calculation. Section VI is summary and discussion. 

II. REVIEW OF SCHWINGER PROPER TIME FORMULATION FOR DIRAC 
OPERATOR WITH A HARD FERMION MASS 

Schwinger proper time formulation 0fj] is a kind of regularization method which can 
covariantly regularize the ultraviolet divergence of Feynman loop diagrams. It is especially 
effective in dealing with the problems related to chiral gauge theories, in which appearance 
of 75 matrix prohibits the naive use of conventional dimensional regularization scheme. The 
most important application of Schwinger proper time formulation is to calculate fermion 
determinant which plays a key role in Feynman loop diagrams from the discussion of last 
section. 

The real part of fermion determinant is related to logarithm of a hermitian operator 
(D^ + m)(D + m) 
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RelnDet(L> + m) = ^Tr \n[(D j + m)(D + m)\ 



(12) 



Take 



E - V 2 = D^D + D f m + mD (13) 

and with help of following relation for a matrix (D^ + m)(D + m) = E — V 2 + m 2 , 

(L> f + m)(L> + m) 



-{D^+m){D+m)T 



i 



-dr = -Ei( 



-7 - ln( 



A 2 ' 
(D^ +m)(D + m) 
A 2 



g (-l)" ( (Dt + m )(p + m ) )n 



n=l 



rain 



A 2 



The real part of fermion determinant in conventional Schwinger proper time formulation is 



Re In Det(-D + m) = - lim 

2 A— >oo 



X) 



(15) 



7 + i n A 2 _ J d * x r *L tr e - mV (x|e^( £ - v2 )| 
where tr is trace for internal symmetry indices, such as spinor, flavor indices, etc. And 

£ = -2ms - 2zm^ 75 + -[7^, 7^7^ + 7 m gP(s - ip7 5 ) + *7 M [0^75(5 - «P7s) + (s - ^75)^75] 
+s 2 +p 2 - [s,p]«7 5 (16) 
7^ = «[V M , V„] = - ^0^75 + Vfw -i{a^ a v ] 
d^f = dpf - i[v^ f] = d^v v - d u v^ - i[v^ v u ] . 

In Schwinger proper time formulation, the regularization is done at lower limit of proper time 
integration r. Since proper time is a parameter irrelevant to symmetry, this regularization 
has the feature of keeping symmetry of the system. 
With help of standard Seely-DeWitt expansion [||], 

1 E 1 



(x\e 



-t(E-V 2 ) 



x) 



1 

To^ 2 



is ~ ~ + ^e 2 - i[V M , [V", E]] - ^-Tl^ - -E* 



6' 



12 



6 



(17) 



+^{E[W, [V„ E]] + [V, [V„ E))E + [W, E] [V M , E)} + T -E' + ■ ■ 



the proper time integration of r in ( [Of) can be finished, substitute ( fL6|) into result formulae, 
we get the expansion of Re lndet(.D + m) with hard fermion mass m, 



RelnDet( J D + m) 
= RelnDet(/9 + m) 



1 



32tt 2 a^oo 



lim / d x tr 



m 



8m[A 2 + m 2 (ln — + 7 - l)}s 



5 



77? 4- 2 777 

_ 8m 2 (ln — + 7 )a 2 - -[Vf - -(In — + T + l)(d^a u - d v a^a v - d"a") 

^ + 7) + y]a 4 + [^(ln ^ + 7) + |] V^a* - 4[A 2 + m 2 (3 In ^ + 3 7 - 1)] S 2 

77?^ '/TO ^ 2 777 ^ 

-4[A 2 + m 2 (ln — + 7 - l)}p 2 + (16 In — + 16 7 + lQ)msa 2 - -(In — + 7 )K^ 
A z A z 6 A z 

Q 777 ^ 1^ 777 ^ 777 ^ 

+i[-(ln ^ + 7) + + 8m(ln — + 7)3**% + e^[-2(\n — + i)V ap d^a u 



m 



+4z(ln — + 7 + 2)a a a p d l) ,a v ] + 0(p 



(18) 



where tr/ is the trace for flavor index and we have isolated out the external fields independent 
term 



RelnDet(/9 + m) 



— lim / d 4 x / — tr e 

2 A^oo 

1 



(x\e Td \x) 



32tt 2 a^oo7i r 3 



00 dr 



lim / — e 



d x trl 



I 777 2 /" 

lim [2A 4 - 4A 2 m 2 - m 4 (2 In — + 2 7 - 1)1 / d 4 xtil (19) 

32tt j a^oo A^ J 



and the expansion terms in fllBj ) are arranged according to its momentum power in which 
6^, Vp, a M are treated as order p and s, p as order p 2 [f|. Hermitian matrix 75 is defined 
as 75 = — 7°7 1 7 2 7 3 which leads relations tr/(7 5 7 At 7^7°7 p ) = — 4e^ CTp and ^^^f- 



T9 



p fia 



There are two ways to calculate coincidence limit of fermion propagator, one is direct 
perform functional differential with respect to corresponding external fields and the other 
is taking similar computation procedure as that for fermion determinant. We will see two 
methods give exactly same results. 

First we consider the functional differential of fermion determinant. Use identity 



<5Tr m(D + m) 

5J u p(x) 
5Tr ln(Dt + m) 



J d*yd 4 z {D + m)-^'°\y,z) 5D ° j l^ ) = (D + m) _1 ' p<J (x, x) 



J d 4 yd 4 z (D^ + m)- 1 ' p ' tT '(y,z 



5D^' p \z,y) 
8Ji a P(x) 



(7?t 



m 



.-l,pcr 



(x,x) (20) 



where a = (a£), p = (b() (a, b are flavor indices, £, ( are Lorentz indices). The scalar part of 
coincidence limit of fermion propagator then is 



5ReTr ln(D + m) 



Ss ab (x) 



5Tt \n(D + m) <5Tr ln(D* + m) 



Ss ab (x) 



5s ab (x) 



d A y 



5Tr \n(D + m) 5J ap (y) 5Ti lnpt + m ) 5J^ p {y) 



SJ a P(y) 5s ab (x) 



6J^P(y) 5s ab (x) 



x, x) + (D' + m) 



-1,(60(00, 



x, x)] 



-. — 9 lim 

47T A— >oo 



2 2 

m[A 2 + m 2 (ln ^ + 7 - l)]5 6a - [A 2 + m 2 (3 In ^ + 3 7 - l)]s ba 



m 



+ (21n^ + 2 7 + 2)m(a 2 ) ba 



+ 0(p 4 ) , 



(21) 



similarly the pseudo scalar part is 



.5ReTr]n(D + m) 



Sp ab (x) 



lim 



<5Tr \n(D + m) 5Tr lnpt + m) 



5p ab (x) 



m 



+ 



5p ab (x) 



m 



47T 2 A->oo 

the vector part is 
.(SReTrlnm + m) 

^ 

<5^(x) 

= o hm 

32tt 2 a^oo 

r 8 „ m 2 



i[A 2 + m 2 (ln — + 7 — l)]p ba + im(ln — + 7)^%) 



+ 0(p 4 ) , (22) 



5Tr ]n(D + m) <5Tr ln(L>t + m ) 



5v ab (x) 



5v ab (x) 



8i 



[d u a u , 
16, 



Pi? 771 ^ $ 771 ^ 

-(In — + 7 + l)[(d"a" - (fa"), a,] - -(In — + 7 )cfV^ 



2 2 

+i[|(ln ^ + 7) + yR(K, a"]) + 8*m(ln ^ + 7 ) [p, o"] + 2ie"«*(ln ^ + 7) <2id»d a a f 



+K,v c 



<yp\ 



+ 0(p 5 ) , 



(23) 



the axial vector part is 
.<5ReTrln(D + m) 



5af(x) 



— r lim 

47T A^oo 



5Tr \n(D + m) 5Tr \n(D* + m) 



1 



5af{x) 
1 m 2 



2m 2 (ln — + 7)0" + - [d^d v a v ] - -(In — + 7 + l)cf (cPa„ - d v a») 



A 2 



A 2 



-[^(ln ^ + 7) + ^](aV + a"a 2 ) + [|(ln ^ + 7) + \}a u a^ 

Tn 2 I ?t? 2 2 

+2(ln — + 7 + lMsa" + a^s) + i[-(ln — + 7) + - V^a,) 

A" 1 d A z d 



777^ I 777^ 

-m(ln — + <i)d»p - -(In — + 7K V„ 



6a 



+ 0(p 5 ) 



(24) 



Since tensor fields are not introduced into the theory, we cannot obtain the tensor part by 
performing functional differential. 

We then take the second method, write propagator in presence of external fields as 



(D + m) 1 (x, x) 

= [(D ] + m)(D + m)]-\D ] + m) 



(x,x) 



7 



lim / dr e- m ' T (x\e- r{E - v2 \D^ + m)\x) 

- (x\e- r{E - v2) f\x} + (x\e- T{E - v2) \x)[m-s(x) - ip(x)j 5 



lim 

A^oo J J, 

A 



A- 



dr e 



-2i<fi (ar)7 5 ] 



Similar as expansion ([T7|), we have 



[x 



16tT 2 



+ 



(25) 



(26) 



Substitute ( p^ ) back to ( p5|) and finish the integration of r, we obtain the expansion of 
(D + m) _1 (x, x) with hard fermion mass m, 

(D + m) _1 (x, x) 

— — lim <^ [A 2 + m 2 (ln — + 7 - l)](m - ijyy 5 ) - 2im 2 (ln— + 

l07T z A^oo L A z A z 



75 



n? 



m 



m 



[A + m 2 (3 In — + 37 - l)]s + z'm(ln — + 7)gZ%7 5 + 2m(ln — + 7 + I) a 



A 2 



A 2 



+ 



7 777 777 1 

-(In — + 7^%, + m(ln — + 7) (pa" - a"p) - -[(eU"K - a»(d u a u )} 



I 777 

+-(ln — + 7 + l)k(dV - cf a") - (ePa" - dVK 
3 A 2 



m 



+ -(ln — + 7 + 2)d u ( afl a u - a,a M ) + e/H-(ln — + 7)^4 



2 2z 



I 771 3 777/ = 

+ 4 (ln A 2 " + 7 ~ 2)a ' Kp " 4 (ln A 2 " + 7 + 2)V(7paV + J a ° a ^ 



m 



OTi(ln — ■ + ^)d^p 
A 2 



2 • 2 

777 % 771 

-2im(In + 7 + iX^s + sa M ) - -(In — + 7 + l)d u {d u a^ - d,a v 
A z 3 A z 



m 



+|(ln ^ + 7 + 2){a v V^ - V^a") - -d^a u + -(In — + 7 + 4)(aV + a"a 2 ) 

2z 777^ % Tfi^ I 1 

--(In — + 7 + 2)fl,aV - e/ 7H--(ln — + 7 KK P - -a a d v a p + -{d a a p )a v } 

+tensor terms + 0(p 4 ) 



7s7 



(27) 



One can easily check that above result can exactly reproduce results (|2l|), (p2|) ,(p3D and 
(0). So two kinds calculation give same results for coincidence limit of fermion propagator. 



III. GENERALIZED SCHWINGER PROPER TIME FORMULATION FOR 
DIRAC OPERATOR WITH DYNAMICAL FERMION SELF ENERGY 

In quantum field theory (|5|) without bare fermion mass, the Lagrangian is chiral symmet- 
ric. Usually the realization of chiral symmetry of the system is characterized by symmetry 
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order parameter (0\tpip\0) . If (Ol^^lO) 7^ 0, we say the chiral symmetry is dynamically bro- 
ken. This happens only when fermion has a nonzero self energy E(/c 2 ) which is in general 
real and momentum dependent. The value of S(£; 2 ) can be determined through solving 
Schwinger-Dyson equation. If we need to use a phenomenological Lagrangian to represent 
this DCSB effect, since momentum dependent fermion self energy E(/c 2 ) in coordinate space 
is represented by £(— d 2 )5(x — y), the naive generalization of fermion interaction for ex- 
plicit chiral symmetry breaking Lagrangian ip{D + m)ip to DCSB Lagrangian should be 
ij)\D + E(— d 2 )]ip. We argue this is not suitable since original ip(D + m)ip is invariant under 
Uv{Nf) local symmetry transformation, 



ij)(x) - 


-> ip'(x) 


= h)(x)ij}(x) 


s{x) - 


-> s'(x) 


= h' (x)s(x)h(x) 


p(x) - 


-> p'(x) 


= h\x)p(x)h(x) 


v^x) - 


-> v'^x) 


= h'(x)v fl (x)h(x) + fv (x)[idph(x)] 


a pix) - 




= h\x)a fl (x)h(x) , 



(28) 



which leads to 

D x + m — ► + m)' = <j) x — if \x) — Zji '(^75 + m — s'(ar) + ip\x)^ 

= W(x)[(jl x — if (x) — i$ (x)75 — s(x) + ip(x)^5 + m]h(x) 

= h\x){D x + m)h(x) (29) 

Here the covariance is due to the property of constant m which leads 

h'(x)mh(x) = m , (30) 

and then ip{D + m)ip is invariant 

^(D + m)i) -> ~$'(D + m)V' = ~$hh){D + m)hh^ = ${D + m)^ (31) 

If we change m to E(— <9 2 ), relation (|2~9"D is no longer valid due to differential operator 
dependence of E, 

/i f (x)E(-9 2 )/i(x) = Ef-Zi^x)^^^)] = E - (fy + ^(^[^(x)]) 2 ^ E(-<9 2 ) . (32) 

So, naive generalization to DCSB Lagrangian ip[D + E(— d 2 )]i/j is not invariant on Uy(Nf) 
symmetry. To implement this local symmetry ()2H), instead of considering E(— d 2 ), we need 
to consider E(— V ) in which 
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V" = d" - iu^x) , (33) 

the bar over V M is to specify the difference of present derivative with that introduced in ([[]). 
Use (p8|) , we find V M transform as 

v: - v:' = ^ - ^(x) = h\x)Xh{x) . (34) 

Then 

E(-Vj) - S(-V^) = E[-^(x)^(x)] = /J(*)£(-V>0r) (35) 

and 

[A* + E(-VS] ^ + E(-V^)] = fct(x)[D a + £(-V 2 )]%) • (36) 

So, ^[£) + £(— V 2 )]ip is invariant under transformation ( p8|) and can be treated as correct 
generalization of t/>(.D + m)ip. Note , in principle, one can add in V M other terms which are 
covariant under transformation (p8|), such as a M 75 multiply by a constant. For these kind 
terms, symmetry itself is not enough to fix them completely, therefore they can be treated 
as some extra interactions which are beyond our choice of V' 1 and £(— V 2 ). Our choice of 
£(— V ) is the minimal generalization required by symmetry to incorporating in fermion self 
energy into the theory. 

The generalized real part of fermion determinant now is 

RelnDetfD + S(-V 2 )] = -Trln + £(-V 2 )][Z> + S(-V 2 )] 



= ii m [°° — Tre -r[B-v 2 +E 2 (-v 2 )+jg(v 2 )+9(v 2 )^-,is(-v 2 )] (37) 

2A^ooJi T V ' 

where 

E - V 2 + S 2 (-V 2 ) + Jg(V 2 ) + g(V 2 )K - 4 S(-V 2 ) = [£>* + S(-V 2 )][,D + S(-V 2 )] (38) 
and 

[4 S(-V 2 )] = 7^[rf M S(-V 2 )] = 7 M (5,S(-V 2 ) - i[v„ S(-V 2 )] 

E = ^[7 M , + 7^(s - ipj 5 ) + i^[a^ 5 (s - ip-f 5 ) + (s - ipTs) Vft] + s 2 + p 2 - [s,p]«75 

#(x) = #(x) = E(-x) J = -ifi 75 - s - zp7 5 if = -zji 75 - s + ipj 5 ■ 

The generalized coincidence limit of fermion propagator is 
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[D + Z(-V 2 )]-\x, 



x 



[ D \ + H(-V 2 )}[D + X(-V 2 )]) [£> f + S(-V 2 )]] (x, x) 
= lim / rfr (x|e-^- v2+s2 (- v ) +J ^ v )+^ v )*H<e(-v )]r_y _ _ ip{ , 

A 7 

-2z^(a;)75 + S(-V 2 )]|x) (39) 
For safety of further calculation, we limit ultraviolet behavior of S(/c 2 ) satisfy constraint 



k 2 







(40) 



For fermion determinant and coincidence limit of propagator, (|37|) and (^) tell us that the 
key now is to calculate matrix element of 



-r[E-V 2 +S 2 (-V 2 )+J 9 (V 2 )+§(V 2 )A'-^S(-V^] 



{x \e 



and 

(x\e 



-r[£-V 2 +E 2 (-V )+Jfl(V )+fl(V )A--jiE(-V 



j;> 



-y — s(x) — ip(x)"f 5 



(41) 

(x) 75 + S(-V 2 )]|x) 



which are much more complex than (|i~7|) and (^). Since now the operator on the exponential 
is beyond original second- order elliptic partial differential operator, £ in principle include 
arbitrary high order of differential operator and except the constraint ([5]), detail function 
S on the differential operator is still unspecified. 

In original Schwinger proper time formulation, there are two methods to calculate matrix 
elements ([H]) and (pE|). One is based on recursion formula and the other is momentum space 
calculation. Recursion formula can be obtained only if we can obtain the matrix element 
with vanishing external fields which now is impossible, since unspecified detail momentum 
dependence of S(fc 2 ) prevent us to obtain it. In this work, we take second method to calculate 
the matrix element in momentum space, 

( x \ e -T\E-V 2 +&{-V 2 )+Jg{V 2 )+~g(V 2 )K-4 E(-V 2 )] ij 

d A k 



(2tt)< 



exp < — t 



E(x) - Vi - 2ik • V x + AT + Ir(-V - 2ik -V x + k' 



+Jg{v x + 2ik ■ V x - k 2 ) + g(y x + 2ik ■ V x - k 2 )K - i S(-V 2 - 2ik -V x + k 2 



• (42) 



Assigning E, J, K to be order of p, we can take low energy expansion for (fl2f). The result is 
given in appendix B. Similarly 



(x\e 



-r[B-V 2 +E 2 (-V )+Jg(V )+g(V )K-4 £(-V )] [ 



-y — s(x) — ip(x)7 5 



(x) 7 5 + S(-V 2 )]|x) 
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d 4 k 



cxp 



E (x) - V 2 X - 2ik ■ V x + k 2 + S 2 (-V - 2ik -V x + k 2 



+Jg{V x + 2ik ■ V x - k 2 ) + g(V x + 2ik ■ V x - k 2 )K - i £(-V; - 2ik -V x + k 2 



x [-f x -ij/t - s(x) - ip{x) lb - 2%i (x) 75 + £(-V* - 2%k ■ V x + k 2 )]\x) 



(43) 



We do not write down the result of expansion for this matrix element explicitly. Instead we 
give the final result in terms of external fields s,p, t^a^ for coincidence limit of fermion 
propagator and discuss it in section V. 

We must note that the difficulty of present calculation is that the momentum integration 
now can not be finished due to unknown behavior of S(A; 2 ). While the chiral covariance of the 
result rely on the achievement of momentum integration in original Schwinger proper time 
formulation. We need to find a way to keep the covariance of chiral symmetry before finishing 
the momentum integration. Fortunately, we found that those non-covariant terms are all 
reduced to total divergence terms at momentum space which with constraint ( fTO|) vanish. 
So with invention of these total divergence terms, we can still obtain a chiral covariant result 
as that in conventional Schwinger proper time method even before we analytically finish the 
momentum integration. 

IV. FERMION DETERMINANT WITH MOMENTUM DEPENDENT FERMION 

SELF ENERGY 

We can parametrize the result fermion determinant as 



RelnDet[ J D + S(-V 2 )] 



d xtrf 



C s + da 2 + C 2 [d^} 2 + C 3 (d»a u - d u a^)(d^a u - d u a^) + C 4 a 4 



(44) 



+C 5 a^a u a,a u + C 6 s 2 + C 7 p 2 + C 8 sa 2 + CgV^V^ + doV^ a„a v + Cupd^a 



+ 0(p 6 



Substitute the definition of E, J, K, g and g into result (|B1| ) given in appendix B, we obtain 
coefficients Ci which are related to E by 



C = -4 / dkZ k X, 



Ci = 2 / dk 



'-2Y? k + k 2 H k Y,' k )X 2 k + (-2E 2 + k 2 ^)^ 
X 2 , Xu k 2 „, X, 



Xk 
A 2 
k 2 



C 2 = 2j dk[2A k Xl + 2A k ^ + A k ^ + ^' k 2 ^ + k -Y! 2 k X 2 k 



C, 



dk 



2B k X k + 2B k 



X\ 
A 2 



Bi 



X k fc 2 „,oX 



A 4 



+ 



A 2 



i v'2 V-2 
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-2 I dk 

'2 



+ + 18 )( * ^ + Iv^ 



X k . 
A6' 



(4E 



-)(2X, 3 



2 A"j„ AY . At, r. 



A 2 



A 4 

-4S 4 fc 



A 2 



2,k Sj, k a 
V L+ 18 )(6Xfc + 



A 2 A 4 A 6 fcv A 2 



'A 4 " 



A 2 " 



C 6 

Cu 



-2 
-2 
-4 



(3E 2 - 2A; 2 E fe E' fc )X 2 + [2E 2 - k 2 (l + 2E fc E' fc )] 



ATfc 
A 2 



(E 2 - 2k 2 Y> k Y>' k )X 2 k - k 2 (l + 2E fc E 



(4^ 



A; 2 E fc )X| + (4E 3 fc + A; 2 E fc )§ + (2E| 



A 2 



1 . <>. . \ Xh 
-k 2 Z k )-i 

2 A 4 



„ Xu 



-3E fc X 2 



-h 2 y' y" 

^k^k 



h 2 y k ^)x k 



Xk 

-Ck + Dk)-^2 ~ ^ k 



D k )X 2 k + 2E k X i 



where 



U / dk 



XI 



4 / dk 



Xk k 2 
A 4 ' 2 
1 



2F k Xl + 2F k ^ + F k ^ + ^' 2 § + ^X? 



k 2 
2 



Vk-\k 2 K)^ + {Y>k--k 2 V k )X k 



dk 



d A k 

y(k 2 ) 



TV 



Xu 



1 



A; 2 + E 2 (A; 2 



(45) 

(46) 
(47) 



and A k , B k , C k , Dk, E k , F k depending on E(fc 2 ) are given in appendix A. We see that asymp- 

fc 2 +S 2 (fc 2 ) 

totic behavior of E(& ) ( f40|) insure factor e ^ appeared in integration measure ( f46|) 
is an ultraviolet damping factor which will keep our momentum integration convergent. 

fl4"4| ) and fl45| ) are our final result for the real part of fermion determinant with presence of 
dynamical quark self energy. The result in this paper is only up to order of p 4 , one can easily 
generalize the calculation to higher orders of the momentum expansion. As a self check of 
theory, take S(/c 2 ) be constant m, in the limit of A 2 — > oo, the momentum integration in 
(fj5|) can be finished, the result gives 



Ci 



E=m N r , , 9 , . , m 2 

— - ~^-2 m l A + m ( in a 2 " + 7 " 1)] 

2=m N c 2 m 2 

— ^4^ m ( in A2 +7) 



13 



c 2 - 

c 4 - 
c 5 - 
<V 
<V 

Cs- 

^10 



s= m AT C m 2 
--^48^ ( '"a? +1 + 1) 



24tt 2 v A 2 



iV r „ m 2 

-Xln—+ 1 + 2) 



24tt 2 v A 2 



8tt 2 



m 



A 2 + m 2 (3ln— + 3 7 - 1) 



A 2 



m 



A + m (In— + 7 - 1) 



s=m JV C m 2 

+ 7 + 1) 

s=m A r c m 2 



48tt 2 v A 2 
iJVU m 2 



12 ^ + " + 2) 

s=m A^ c m 2 
Cu >--^m(7n— +7) • 



(48) 



Substitute them back into (0), we reproduce original result fll~8|). So our generalized for- 
mulation can easily recover the conventional Schwinger proper time result. 

V. COINCIDENCE LIMIT OF FERMION PROPAGATOR WITH MOMENTUM 

DEPENDENT FERMION SELF ENERGY 



The computation of fl39|) shows that the scalar, pseudoscalar and axial vector parts of 
coincidence limit of fermion propagator are 
1 



-(1)« 
2 V ; 



r2y|-l,(6C)(a€)| 
6a 



[D + ^(~^)]- 1 ' {Km) (x,x) + [rf + S(-V T )] 



12x1-1,(60(00. 



x, x) 



C + 2C e s + C 8 a^ 



+ 0(p 4 ) , 



2 (7s) ^ 



[D + E(-V )]- 1,(6C)(a0 (a;,a;) - [L> f + S(-V T )] 



t2xi-l,(6O(a0, 



X, x) 



2C 7 p + Cud^a. 



ba 



+ 0(p 4 



2(757m) CC 



[D + £(-V / )]- 1 »^)^)(a;, x) + [£>t + S(-V T )] 



t2x 1 -l,(6O(a0, 



X, X) 



(49) 



(50) 



2C 1 a /i - 2C 2 \<Fd v a v \ + 2C 3 d u (d»a u - d u a») + 2C 4 (aV + a^a 2 ) + AC^aTa" 



+C 8 (sa" + a»s) + C 10 (a^^ - V^a v ) - C u d"p 



ba 



+ 0(p 5 ) , 



(51) 
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where the coefficients Ci are given in last section. Note the computation is done for fermion 
determinant and propagator independently, and for scalar, pseudoscalar and axial vector 
parts, we recover the result of which we directly performing functional differential with 
scalar, pseudo scalar and axial vector fields for determinant. Unfortunately, this property 
no longer valid for vector part of propagator. We parametrize the corresponding coincidence 
limit of fermion propagator as 



2(7, 



[D + S(-V 2 )]- 1 ' (6C)(< ^ ) (^^) " \D ] + S(-V t2 )]- 1 ' (6c)( ^ ) (^^) 



ba 



2iC 2 [d v a v , a"] + 4iC 3 [d^a u - d u a" } a v \ + 4C 9 d u V^ + C w d u [a^ a v \ + iC u \p, a A 
+0(p 5 ) , (52) 

The coefficients d with overline on them are to specify the difference with their original 
coefficient Ci from functional differential with vector fields. We list down their values in 
appendix C 

Now we discuss the reason to cause this difference. With momentum dependent fermion 
self energy, the coincidence limit of fermion propagator can not be obtained directly as that 
in hard mass case from the functional differential of external fields for fermion determinant. 
Since the formula corresponding to (p0|) now is changed to 

<5Trln[£) + S(-V 2 )] 
8J a P{x) 

= {D+ E(-V 2 r'-(x,x) + / d'yd'z \D + n-^)\- 1 ^b^) m ~j)^ Z ' y) (53) 

and its hermitian conjugate formula. The second term in ([53]) is a new term which rely on the 
momentum dependent S(A; 2 ) and vanishes when self energy is independent of external fields. 
Because of this extra term, functional differential of external fields for fermion determinant 
no longer give propagator. Note that £(— V ) in the second term is independent of scalar, 
pseudo scalar and axial vector fields and the corresponding differential with scalar, pseudo 
scalar and axial vector fields vanish, so this extra term donot contribute to scalar, pseudo 
scalar and axial vector parts of coincidence limit of fermion propagator. This explains why 
direct calculation of scalar, pseudo scalar and axial vector parts of coincidence limit of 
propagator can recover result of functional differential for determinant. 
Paramitrize the second term as 



15 



[D +h{ V )\ (y,z) sjmm){x) 

2tC 2 [d u a u , a"} + AiC 2 \d t "a u - d u a», a u ] + AC 9 d u V^ u + C lo d v [a", a u ] + iC u [p, a"] 
+0(p 5 ) , 



ba 



(54) 



We can use our generalized Schwinger proper time method directly calculate these coeffi- 
cients. 



Tr 



|fl + E( _vY^ + E <- v2 >l 



5J a P(x) 



Tr 



-i 



[£t + s( _v 2 )][D + s(-v 2 )]) >t + S (-v 2 )]^ + s( - v )] 



8J°p{x) 



lim / dr Id y tx{y\e 



-r[£-V 2 +E 2 (-V )+Jg(V )+9(V )K-fE(-V 



[-f - s(x) - ip(x)j 5 



(55) 



8J a P{x) 

Combine above result and its hermitian conjugate together, we can calculate coefficients C 
in (|54|). We list down the result coefficients C in appendix C. One can check that summation 
over C and C coefficients together recover original C coefficients: 



Qf A ~\~ Cj 



i = 2,3,9,10,11 



(56) 



VI. SUMMARY AND DISCUSSION 

In this paper, we have generalized conventional Schwinger proper time method for stan- 
dard Dirac operator to incorporate dynamical chiral symmetry breaking. The physical 
output of this generalization is the real part of logarithm of fermion determinant and coinci- 
dence limit of fermion propagator in presence of momentum dependent fermion self energy 
X(/c 2 ). The mathematical progress of this generalization is that the operator on the expo- 
nential in the key matrix element of Schwinger proper time formulation is now generalized 
from original second-order elliptic differential operator E — V 2 (see (|i~7|)) to arbitrary high 
order differential operator E - V 2 + £ 2 (-V 2 ) + Jg(V 2 ) + g(V 2 )K - i £(-V 2 ) (see ©). 
This high order differential operator dependence is represented by an unspecified function 
£(— V ) which physically is fermion self energy and characterize the dynamical chiral sym- 
metry breaking of the system. The generalized formulation automatically keep the local 
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symmetry Uy{Nf) given by fl28|) which is residual symmetry after explicitly breaking of 
local chiral symmetry U^Nf) ® Un(Nf) by constant fermion mass. 

As discussed in the introduction of this paper, one can further generalize this local 
symmetry Uv(Nf) invariance of the formulation to the invariance of full chiral symmetry 
UL{Nf)^)Uji(Nf). The price is that the symmetry must be realized nonlinearly, i.e., we need 
to introduce into theory a local field Q(x) which transform nonlinearly under local chiral 
symmetry (|3|) as (|^) and rotated external fields defined as (^|). Replacing the external fields 
in D + S(-V 2 ) in © and (||) with D Q + S(-V^), = - iv&(x). One can easily 
check that under local chiral symmetry Ui(Nf) <g> Un(Nf) transformation (||) and (0), we 
have (0), (0) and 

D n + E(-Vn) ^D' n + S(-Vn) = h\x)[D n + E(-Vl)]h(x) (57) 

i.e., they are covariant quantities. In terms of rotated quantities, once our formulation keep 
this local hidden symmetry Uy(Nf), it automatically keep original local chiral symmetry 
U L (N f )®U R (N f ). 

In this paper, we do not calculate the imaginary part of fermion determinant and fermion 
propagator are only limited to its coincidence limit, its general nonlocal part is not com- 
puted. These are all under investigation and we will present results in separated papers. 
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APPENDIX A: COEFFICIENTS DEFINITIONS 

A k = ^ 2 E fc E' fc (-l - 2E fc E' fc ) - 1e 2 (-1 - 2E fc E^) + ^ 2 E 2 (E', 2 + S fe E' fe ') 

+^ 4 (e' 2 + s fe E£) 

B k = h^^i-l - 2S fc S' fc ) - l -Y? k {-\ - 2E fe E' fc ) + ^ 2 E 2 (E', 2 + S fe E' fe ') 

+^ 4 (si 2 + s fc s' fc ') + h\-i - 2^' k ) 

Ck = g - gS fc S' fc + -k 2 T! k 2 

D k = -\k 2 V k 2 + h^i-l - 2E fe E' fc ) + ^ 4 E' fe E£(l + 2E fe E' fc )] - h^i-H' 2 - E fc EZ) 

+ -£; 2 E A .E / fc (-E / fc 2 - E fc E' fc ') 
E k = ^ 2 E fe E' fe (-l - 2E fe E' fc ) 2 - h^il + 2E fc E^) 2 

^ = \k 2 ^' k - \k 2 {^' k ) 2 - h 2 + ^E' fc + ^ 2 E 2 (E' 2 + E fe E») + ^ 4 (E' fc 2 + E fe E») 
+ ^ 2 (-l-2E fe E^) + lA; 2 

APPENDIX B: GENERALIZED SEELY-DEWITT EXPANSION 

The detail calculation gives: 

( x \ e -T[E-VHsH-V 2 )+Jg(V 2 )+~g(V 2 )K-4 £(-V 2 )] |^ 

= / ^' Tf {^ + ^ + ^ + \k 2 {- T -f + r 2 ff" - y/' 3 )[V^, Vl[V„ V,] 
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-^f 'A: 4 75([V M , V]]a„ + a^, [V„ V]]) + i^f f"k%(-2[a\ [V\ [V M , V,]]] 
+ [V M , V^aJ + [V„ajV M , V] + [V„, [V M , V>,] + [V„ a,[V\V v }]) 

-^/' 3 75([V„a,][V M , Vl + [V^, Vl[V„a,]) - V[V M , [V„a 2 ]]) 
-^/"A; 4 ([V M , a M ] [V u , a v \ + 2[V", o"] [V M , a,] - 2[V M , [V", a,a v \\ - 2a,a u [V\ V"]) 
+ ^f 2 k 4 (a»[V^ [V M ,a„]] + [V M , [V^K + [V M , a"] [V M , a,] + [V", a"]^, a„] 

-a^[V M , V"] - [V M , V> M a„) - ^A 2 75 (K [V M , a 2 ]] + [V M , [a M ,a 2 ]]) 

?V 4 _ _ _ 

+ 3gAS(^,a1k,flM] + k,a,][V^,a1 + [V M , [a„a 2 ]] + K [V M ,a 2 ]]) 

-J/V^, [V M , F]] + ^/VfV", [V M , F']] + ^A; 2 (/ J[V*, [V„ J]] 

+^ 2 [V M , [V„ K]]K + g'g'[V", [V„ FJ]]) - A; 2 (2F'[V M , [V M , F]] 
+2[V M , [V M , F']]F + 2[V M , F'] [V„ F] + 2F[V M , [V M , F']] + 2[V M , [V„ F]]F' 
+2[V M , F] [V„ F'] - [V M , [V„ J<?'F + F£'F]] - [V M , [V„ F]] J«/ - ~g'K[V^ [V„ F]) 

+^A 2 (F[V A \ [V„ F]] + [V M , [V„ F]]F + [V", F] [V„ F]) 

+ ^k 2 l5 [g'(a,[V> 1 , J] + [V^, a^J] - J[V M , aj) + ^'([V M , a M ]F - [V M , Fa M ] - [V M , F]a M )] 

A 2 7 5 ([V M , [a„ F]] + [a„, [V M , F]]) - i 2 75 (F'[V^, [a„, F]] + [a„ [V M , F']]F 
+[a„, F'] [V M , F] + F'[V M , [a M , F]] + [V M , [a„, F']]F + [V M , F'] [a„, F]) 

"12 _ _ 

+ [V M , Fa M J] + S F[V M , J] + [V M , FJ]^) + [^(-^[V M , FF] - [V M , F]Fa M - [V M , Fa^F] 

-a M [V M , F]F - a„F[V M , F] + [V M , Fa M ]F - F^[V M , F] + [V M , F]Fa M )] 

+ ^A 2 7 5 (F[V M , [a„F]] + [V M , [a,75,F]]F+ [V M ,F][a„F] +F[a„ [V M ,F]] 

+ [a„, [V M ,F]]F+ [a M ,F][V M ,F]) - £ £'A: 2 (F'[V M , H?, V M ]] - [V M Hy , V"]]F') 

^^S'A; 2 ((^[V^F][y,V^+^[y,V M ]][J,Vl) + yS'A; 2 [y,V At ] 2 

+^m 2 ([v,, [y , v m ]f] + [v„F[y, vl] + [[y, vj, [v m ,f]) 

-^S'A; 2 75 ([y,Vla M + a M [y,Vl) + ^S' 75 A; 2 (^[y ! V M ]F + ^F[y,Vl 
+ [y,V^]Fa^ + [y,V M ]a^F + Fa^[y,V M ] + F[y,v M K)j 

+total derivative terms (Bl) 
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--T^9'k 2 lM[^^ J W F ~ «^[V M , F] - J[V", a M F] + [V M , J]Fa M + F[V M , J]a M 



where the total derivative terms are 



total derivative terms 



= 1 ^wA e -% k2r + V - k w + k k2 f 3 + V 2 ^ 

+ ^ 2 75 e-^f [V VV M + V V) + + V%) V 2 ] 

-^A;V^r [V M VV + V 2 V% + a^V M V 2 + a„vV] 

+^£:V r/ / ,2 [V + V V) + KV M + V%) V 2 + 2V A VV„ + V%) V„ 

+2V M V 2 + 2a„vV] + ^A; 2 e-^(- fc2 )/'(-A; 2 )[V 2 a 2 + a 2 V 2 ] 

8 

+ !I^ e -/ [a 2 (a Mv /i + y%) + KV M + V M a> 2 ] - ^e^^a, + a^) 2 

o o 

-^ 2 e- r/ /'[(V% + a^) 2 + 2V"(V"a, i + a^V M K + 2a"(V'a A1 + a^V^V, + 2a 2 V 2 
+2V a 2 + 2/V aJ - — e - r/ [a 2 (a^V M + V%) + + V"a> 2 + a>"V, + a%K] 

oD 

2 

-^e- r/ (^JV 2 + gV 2 K) + ^-e- T/ f (FV 2 + V 2 F) 
4 8 

+ ^e- Tf [g'(JV 2 F + FJV 2 ) + ^(v'i^F + FV 2 K)} - ^ e - r/ (vV + F 2 V 2 + FV 2 F) 
8 24 

+^^([«% + V M a^]F + F[a"V„ + VXD 

• 3 

-^ e - T/ ([a% + V M a^]F 2 + F VV M + VX1 + F[a^V M + V./JF) 

+ ? V^s'[y,v 2 ] -^^e-^s'([y,v 2 ]F + F[y,v^])} (B2) 

and 



<% 2 J ~ d(k 2 ) 2 J ~ d(k 2 f 
F = E(x) + J(x)g(-k 2 ) + g(-k 2 )K(x) 

F' = J(x)g\-k 2 ) + ~g\-k 2 )K(x) 

2 3 

Y 1 (x) = 1 - r[a 2 (x) + F(x)} + — [a 2 (x) + F{x)f - — \a 2 {x)F 2 {x) + F 2 (x)a 2 (x) 

2 6 

+F(a;)a 2 (x)F(x) + F 3 (x)] + ^F\x) 
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Y 2 (x) = —a 2 (x) \2a A (x) + a^(x)a 2 (x)a^(x) + F(x)a 2 (x) + a 2 (x)F(x) + a^F(x)a^(x)} 

-^[F 2 (x)a 2 (x) + a 2 (:r)F 2 (:r) + a M (a;)F 2 (a;)a M (a;) + F(x)a 2 (x)F(x) 
+F{x)a' x {x)F{x)a fl (x) + a>"(x)F(x)a^x)F(x)] 

Y 3 (x) = -^[a 4 (x) + a M (x)a v (x)a fi (x)a v (x) + a M (x)a 2 (x)a M (x)] 

APPENDIX C: C AND C COEFFICIENTS 

C coefficients are defined as 

C 2 = 2jdkA k [2X! + 2^ + ^] 
C 3 = jdk [B k (2Xl + 2§ + §) + X 2 k + §] 
C 9 = - j^Cy|f + X 2 ] 
C 10 = 4tjd~k F k [2X 3 k +2^ + ^} 

C u = Aj dk^ k [^+X 2 k ] (CI) 

with 

A fe = i[2E 2 -fc 2 E fe E' fe ] 
B k = l -[k 2 + Y? k + k 2 ^' k ] 
C k = \\l-^ k + \k 2 V 2 ] 

F k = i[-4E 2 + 5A; 2 E fc E' fc + k 2 } (C2) 
C coefficients are denned as 

6 = 2/ <£ [i,(2X fc 3 + 2§ + §) + y S' fc 2 (^f + X 2 )] 
C 3 = Jdk [B k {2Xl + 2§ + ff) - (1 - \k 2 K 2 )(§ + X 2 )] 
C 9 = y dfc [i(fc 2 E' fe - E fe )E' fe 'X fc + £ fe (^f + X 2 k ) + £ fc (2X fe 3 + 2^| + ^)] 
C10 = 4< / dfc [F fc (2X| + 2§ + + y S' fc 2 (§ + X 2 )] 

C n = -2 f dk k 2 Z> k (^ + X 2 ) (C3) 
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with 

A k = l -{-2k 2 H' k + S fc )E fe (l + 2E fe E' fc ) + ^ 2 (2E 2 + P)(E' 2 + EE") - ^E 2 + ^ 2 E fe E', 
5, = ^E fe (-2£; 2 E' fe + E fc )(l + 2E fe E' fe ) + ^(6E 2 + A; 2 )(E' 2 + 

-y (1 + 2E fc E' fe ) - i(A; 2 + E 2 + A; 2 E fc E' fc )] (C4) 
D k = -^ 2 E' 2 + A; 2 (-iE fe E^ + ^ 2 E' fc ED(l + 2E fe E^) + A; 2 E' fc (^ 2 E' fc - is fc )(E' fc 2 + E fc E' fe )] 
^ = A: 2 S' fc (is fc -iA; 2 E' fc )(l + 2E fe E' fc ) 2 

h = S fc E' fe (^ 2 - ^%E' fc + ^E 2 ) + A; 2 ^ 2 + i*; 2 )^ 2 + S fe E^) - ^(1 + 2E fe E' fe ) + \h>] 



22 



